Electromagnetic waves in a Josephson junction in a thin film 
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G\ '. Abstract 



We consider a one-dimensional Josephson junction in a superconducting film 



P^ ' with the thickness that is much less than the London penetration depth. We 

o' 

J-^ ' treat an electromagnetic wave propagating along this tunnel contact. We show 

L^ that the electrodynamics of a Josephson junction in a thin film is nonlocal 

C^ ' if the wave length is less than the Pearl penetration depth. We find the 

integro-differential equation determining the phase difference between the two 

superconductors forming the tunnel contact. We use this equation to calculate 

the dispersion relation for an electromagnetic wave propagating along the 

Josephson junction. We find that the frequency of this wave is proportional 

to the square root of the wave vector if the wave length is less than the Pearl 

penetration depth. 
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I. INTRODUCTION 

The electromagnetic properties of tunnel Josephson junctions are a subject of intensive 
studies over the past three decadestl. A considerable attention is, in particular, attracted 
to the investigation of the SIS-type Josephson contacts. In this case the tunnel junction is 
formed by a thin layer of an insulator. This dielectric layer between two superconducting 
plates can be treated as a transmission line or a parallel plate resonator when the electromag- 
netic properties are concerned. It follows from this approach that an electromagnetic wave 
with a specific dispersion relation may propagate along the SIS-type Josephson junctionH. 

The existence of this Swihart electromagnetic wave results, in particular, in the self- 
induced resonances, usually referred to as Fiske stepsEi. These specific resonances are ob- 
served as peaks in the current-voltage curves of the tunnel Josephson junctions. The study 
of Fiske steps is one of the methods to treat the electromagnetic properties of the SIS-type 
Josephson contacts. Recently, it was successfully applied to investigate the electromagnetic 
properties of the grain boundaries in YBa2Cu307_5 high-temperature superconductorsa. 

The self-induced Fiske resonances are observed in the presence of an external magnetic 
field, Hf.. In this case the phase difference between the two superconductors forming the 
tunnel contact, </?, is a sum of two terms, i.e., y? = (pu + ^v- The first term, ipu., describes 
the increase of (p along the Josephson junction (y-axis). The value of (^h is proportional 
to the magnetic field Hf. and is given by the formula ifH = ~4:TTXHey/^Q, where A is the 
London penetration depth and $o is the flux quantum. The second term, (pv, represents 
the increase of if in time. The value of ipv is proportional to the voltage, V, applied to the 
contact, and is given by the formula 2eVt/h. As a result the tunnel current density 

J =jcsmip = j^smi^-^t yj (1) 

is described by a wave propagating along the Josephson junctional. The amplitude of this 
wave is given by the Josephson critical current density, jc- 

The self-induced Fiske resonances arise when the frequency of the tunnel current den- 
sity wave, 2eV/h, becomes equal to the frequency of the Swihart electromagnetic wave'^'* 



Thus, the dependence of the frequency, cj, on the wave vector, k, for an electromagnetic 
wave propagating along the Josephson junction determines the positions of the peaks in the 
current- volt age curve of a tunnel contact. 

Usually, the dispersion relation uj{k) is determined by the sin-Gordon equation that leads 

tS 



uj = uj, ^1 + fc^A^, (2) 

where 



is the Josephson frequency, C is the specific capacitance of the tunnel junction, and 



is the Josephson penetration depth. 

The Swihart electromagnetic wave corresponds to the limiting case 1 <^ k\j. It follows 
than from Eq. (^ that uj ~ ujjXjk, i.e., this wave is propagating along the Josephson 
junction with a constant velocity 

The linear dispersion relation uo = Csk results in an equidistant set of peaks in the 
current- volt age curves of the SIS-type Josephson tunnel junctiona§Q. 

We can determine the phase difference if{y,t) in the mainframe of the local Josephson 
electrodynamics, i.e., by the sin-Gordon equation, as long as k\ <C la. It means, in particu- 
lar, that for an electromagnetic wave propagating along the Josephson junction with A <^ Aj 
the dispersion relation cu = Cgk is valid in the region Aj^ <^ k <^ A^^. 

Let us now discuss the general case, i.e., the case when the restrictions on the wave 
vector k are given by the inequalities kdo ^ 1 and fc^ ^ 1, where do is the thickness of 
the insulating barrier and ^ is the coherence length. We treat here an SIS-type Josephson 



junction formed by two superconducting plates. In this case the space distribution of </) is 
one-dimensional and the relation between the phase difference </?(?/, t) and the magnetic field 
in the superconductors is nonlocal if 1 <C kX. As a result the function (p{y,t) is determined 
by an integro-differential equationQB, i.e., the electrodynamics of a Josephson junction is 
nonlocal as far as the region of wave vectors 1 <^ /eA is concerned. Using this equation it 
was shownl3ll3 that the dispersion relation for an electromagnetic wave with I <^ k\ takes 
the form 

uj = c, y-. (6) 

The phase velocity of this electromagnetic wave is inversely proportional to the square 
root of the wave vector k. It results in a non-equidistant set of the values of the voltage V 
corresponding to the self-induced resonances in the current-voltage curve.Q 

This effect is most pronounced when considering an SIS-type Josephson junction in a 
thin superconducting film with the thickness d <^ X. In this case the space scale of magnetic 
field variation in the superconductors forming the contact is given by the Pearl penetration 
depthS 

Aeff = ^ > A. (7) 

Thus, the electrodynamics of a Josephson junction in a thin film is nonlocal if 1 <^ kXcs. 
This region of wave vectors is much wider than the one given by the inequality 1 <^ A;A. 

In this paper we consider an infinite one-dimensional Josephson tunnel junction in a 
superconducting film with the thickness (i <^ A. We show that the electrodynamics of a 
Josephson contact is nonlocal if the space scale of variation of ip is less than the Pearl 
penetration depth, i.e., if 1 -C kXes. We derive the integro-differential equation determining 
the phase difference ip{y,t). We use this equation to calculate the dispersion relation uj{k) 
for an electromagnetic wave propagating along an SIS-type Josephson junction. 



^We present the theory of the self-induced Fiske resonances in the mainframe of the nonlocal 
Josephson electrodynamics elsewhere. 



The paper is organized in the following way. In Sec. II, we consider the electrodynamics of 
a long one-dimensional Josephson junction in a superconducting film with the thickness that 
is much less than the London penetration depth. We treat the general case of an arbitrary 
relation between the Josephson penetration depth and the effective Pearl penetration depth. 
We derive the integro-differential equation determining the phase difference distribution 
along the Josephson junction. In Sec. Ill, we apply this equation to calculate the dispersion 
relation for an electromagnetic wave propagating along the Josephson junction. In Sec. IV, 
we summarize the overall conclusions. 

II. BASIC EQUATIONS 

Let us consider a thin superconducting film (xy-plane) with an SIS-type Josephson junc- 
tion along the y-a.xis as it is shown in Fig. |1|. We treat here the case when A ^ ^, in which 
the London equations govern the fields and currents inside the superconductor. Thus, out- 
side the tunnel contact the relation between the current density j and the magnetic field h 
is given byt3 

, 47rA2 

h + rotj = 0. (8) 

Introducing the vector potential A and combining Eq. (1) with the equation 

h = rot A (9) 

we express the current density j in the form 

J = J^(S-A). (10) 

where outside the tunnel junction the vector field S is given by the formula 

S = ^V9 (11) 

and 6 is the phase of the order parameter. 



The quantities j, A and S are nearly independent on the z coordinate in the hmiting case 
of a thin film, z.e., for rf ^ A. Therefore, in order to find the fields and currents we replace 
the superconducting film with the thickness d <^ A by an infinitely thin current-carrying 
sheet in the plane z = 0. The current density j in this plane is determined then by the 
averaging of Eq. (|r^) over the thickness ctllGJ, which results in 

J = 7-^(S-A)<5(z). (12) 

Let us now choose the London gauge, ie., let us assume that div A = 0. Then, substi- 
tuting Eq. (|12D into the Maxwell equation 

47r 
roth = — j (13) 

we find the equation describing the vector potential A in the form 

~^K^\-lK6{z) = \-2^6{z). (14) 

The vector field S is related to the phase difference 

<^(l/) = ^(+0,2/)-^(-0,2/). (15) 

This relation is given by the equation 

rotS = ^^'(y)%)z, (16) 

following from Eq. (0) and taking into account the singularity of the function 6{x,y) at 
X = 0. The vector z is here for the unit vector along the z-axis. 

Applying the continuity equation divj = to Eq. (|12|) we find that divS = 0. Thus, we 
can present the vector field S as a curl of a certain vector field F, namely, 

S(p)=rotF, (17) 

where p = {x, y) and 

F = F{p) z. (18) 
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Substituting Eq. ([T8|) into Eq. (0) we find the equation describing the function F{p) in the 
form 

^F = -^^\y)5{x). (19) 

Note that the scalar function F{p) determines both components of the vector field S(p) 
reducing by this way the complexity of the problem. 

The current density across the Josephson junction ja;(0, y) is a sum of two terms, namely, 
the tunnel and the displacement current densities 

j,(0,y) = [j, siny. + — -Q^\d5{z). (20) 

At the same time it follows from Eq. (|1^) that the current density jxiS^, y) can be written as 

Jx.(0, y) = -^ [SM y) - ^.(0, y. 0)] 5{z). (21) 

Equating the expressions for the quantity jx{0,y) given by Eqs. ( pO]) and (^) we find that 

J, siny; + — ^ = ^ [Sx{0,y) - Ax{0,y,0)] (22) 

Thus, to derive the closed form of the equation describing the phase difference ip{y,t) it 
is necessary to find the functional relation between 

Ax{y) = Sx{0,y)~Ax{0,y,0) (23) 

and if{y,t). We use here Fourier transformation in order to do it and defining the Fourier 
transforms for A(r) and S(p) as 

/d odp 
Aq,p exp{iqp + ipz) — -y (24) 

and 

S(p) = /S,exp(.qp)^. (25) 

Using Aq p and Sq we can present the value of A^. by the integral 



OO TV 



A, = ^ JqdqJd^iS^^- AD expitqysin^), (26) 



where ^ is the polar angle in the {qx, gjy)-plane and 



A, = / A,,, g. (27) 



The next step is to apply Fourier transformation to Eq. (]1^) , which results in 

iq' + p") Aq,, + \-^ Aq = A;^^ Sq. (28) 

It follows from Eq. (^) that the relation between Aq and Sq has the form 

Aq = ^3-— (29) 

^ l + 2gAeff ^ ' 



and thus 



2gA 



eflF 



Sl-Al = — ^-^^ 5^. (30) 

'^ "^ l + 2gAeff '^ ^ ' 

To calculate the Fourier transform 5"^ we take the derivative of Eq. (|T^ with respect to 
y and substitute S^ instead of dF/dy. As a result it comes out that 

ASx = -^v"{y)5{x). (31) 

It follows from Eq. (^) that the Fourier transform S^ is given by the formula 

'5q = ^^ J dy if" {y)exp{-iqy sin i9) (32) 

— OO 

Combining now Eqs. (|32D , (^), (p6D, and (0) we find the integro-differential equation 



describing the phase difference <f{y,t) in the form 

1 av 



io^j dt^ 



where 



+ sin^ = lj J dy'K{^L_]L) ^"(2/'), (33) 



OO 



K{u) = - [ -^^ dv, (34) 

71 J V+\u\ ^ ^ 

' ' 



Jo(f ) is the zero-order Bessel function, and 



C$r 



Note that Eq. ( |33D can be rewritten asB 



(35) 



1 aV , . h f dy' dip 

uij ot^ n J y' — y oy' 



(36) 



in the hmiting case when the characteristic space scale of the phase difference variation is 
much less than X^q. 

III. DISPERSION RELATION 



Let us now consider a small amplitude electromagnetic wave propagating along the 
Josephson junction. The corresponding solution of Eq. (Rl) then reads 



if = ifQ exp{iky — iut), |y9o| -C 1. 



Substitution of Eq. ( ^7|) into Eq. (^) results in the following dispersion relation 



(37) 



u = ujj ^1 + 2k^X,sljlC{2kX,s) 
where 

oo 

/C(x) = / K (u) exp{ixu) du. 

— oo 

The function }C{x) has the following explicit form 



(38) 



(39) 



}Cix) 



1 , 1 + VI - a;2 
= In 



vrVl — x'^ 1 — VT—x^ 



1 arctan . 

TT y/x"^ — 1 



if a; < 1; 



if a; > 1. 



(40) 



Using Eqs. (^) and (^OD we find, in particular, the dispersion relation uj{k) in the limiting 
cases kXes ^ 1 and kXes ^ 1 
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UJ 



ujJl ^^ln(A;Aeflf), if /cAcfr < 1; 

(41) 



^ ujj\/l + klj, if kXcs > 1. 

Thus, for kXcs ^ 1 and klj ^ 1 the frequency of an electromagnetic wave propagating 
along the SIS-type tunnel Josephson junction in a thin film is proportional to the square 
root of the wave vector. 

The dispersion relation cu oc y/k leads, in particular, to a non-equidistant set of the 
self-induced Fiske resonances in the current-voltage curves for the voltages V > Vc, where 

Note, that for a thin superconducting film the relation kXcs ~ 1 corresponds to a wave 
length that is X/d times bigger than the London penetration depth. 

IV. SUMMARY 

To summarize, we have found the integro-differential equation describing the phase dif- 
ference in case of the SIS-type tunnel Josephson junction in a thin superconducting film. 
We apply this equation to calculate the dispersion relation for an electromagnetic wave 
propagating along the Josephson contact. We have shown that if the wave length is small 
compared with the Pearl penetration depth Aefr the frequency is proportional to the square 
root of the wave vector. 
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FIGURES 
FIG. 1. A thin superconducting film with a Josephson junction (thick fine). 
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